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The energy spectrum for an ultracold rotating Bose gas in a harmonic trap is calculated exactly 
for small systems, allowing the atoms to occupy several Landau levels. Two vortex-like states and 
two strongly correlated states (the Pfafnan and Laughlin) are considered in detail. In particular, 
their critical rotation frequencies and energy gaps are determined as a function of particle number, 
interaction strength, and the number of Landau levels occupied (up to three). For the vortex- 
like states, the lowest-Landau-level (LLL) approximation is justified only if the interaction strength 
decreases with the number of particles; nevertheless, the constant of proportionality increases rapidly 
with the angular momentum per particle. For the strongly correlated states, however, the interaction 
strength can increase with particle number without violating the LLL condition. The results suggest 
that in large systems, the Pfaffian and Laughlin states might be stabilized at rotation frequencies 
below the centrifugal limit for sufficiently large interaction strengths, with energy gaps a significant 
fraction of the trap energy. 
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I. INTRODUCTION 

Over the past few years, there has been much inter- 
est in the possibility of coaxing ultracold atomic gases 
into the neutral analog of a quantum Hall state. Various 
approaches have been proposed, which include rapidly 
rotating the confinement potential [1-14], and generat- 
ing effective magnetic fields for particles confined in op- 
tical lattice potentials [15-18] and harmonic oscillator 
traps [19, 20]. 

One of the many intriguing prospects for these systems 
is the possibility of producing the bosonic analog of the 
Pfaffian [2, 7, 10, 14], also known as the Moore-Read 
state [21], which is believed to give rise to the v = 5/2 
fractional quantum Hall plateau in two-dimensional elec- 
tron gases [22] . The quasiparticlc excitations of this state 
are believed to be non-Abelian anyons obeying fractional 
statistics [23-25]; in principle, non-Abelian anyons could 
be physically wound around one another to operate an 
intrinsically fault-tolerant quantum computer, protected 
from errors by the topological properties of the underly- 
ing state [23, 26, 27]. 

One of the fundamental sources of error in these topo- 
logical approaches to computing is the thermal excitation 
of unwanted quasiparticles that contribute to the wind- 
ing: the error rate is intrinsically related to the size of the 
energy gap separating the ground state from the next ex- 
cited one. One might expect that strengthening the par- 
ticle interactions would monotonically increase the split- 
ting between sub-bands in the Landau levels; this would 
also allow the system to access a larger number of quan- 
tum Hall states. Indeed, for ultracold atomic systems 
the size (and sign) of the scattering length (character- 
izing the interaction strength) can be varied arbitrarily 
through the use of Feshbach resonances [28]. It is not 
known, however, how large interaction strengths affect 
the quantum Hall states and the excitation gaps, except 
right on resonance [10]. 

The difficulty is that calculations used to locate the 
various quantum Hall ground states usually rely on the 
Lowest Landau Level (LLL) approximation, which is 



valid only in the weakly interacting regime. This is turn 
requires that the particle densities are low, the interac- 
tion strength is small, or the rotation frequencies are ex- 
ceedingly close to the harmonic oscillator confinement 
frequency. But robust quantum Hall states with large 
excitation gaps require strong interactions. 

In the present work, we exactly calculate the ground 
states for rotating bosons for small systems without em- 
ploying the LLL approximation. The central goals are (1) 
to determine an upper bound to the interaction strength 
consistent with the LLL approximation; (2) to obtain 
the physical parameters (number of bosons, scattering 
length, rotation frequency) to unambiguously select a 
given quantum Hall state; and (3) to show that nega- 
tive (attractive) interactions always lead to instability 
for large densities, contrary to prior claims [29, 30]. The 
main result of this paper is that the widely accepted cri- 
terion for the LLL approximation, which requires that 
the inter-particle interaction energy be smaller than the 
Landau level spacing, is generally too restrictive: we ex- 
pect that the most interesting quantum Hall states can 
be obtained for relatively large particle numbers and ro- 
tation frequencies within reach of current experiments. 

In Sec. II, we describe the Bose gas and our approach 
to calculating the boson energy spectrum. Limits on 
the validity of the LLL approximation are investigated 
in Sec. Ill by considering the effects of higher Landau 
levels on the boson energy spectrum. The case of attrac- 
tive Bose gases is discussed in Sec. IV, and concluding 
remarks are provided in Sec. V. 



II. NUMERICAL METHOD 

The system under study consists of ultracold interact- 
ing bosons subject to a cylindrically symmetric harmonic 
trapping potential, which is rotated around the z-axis 
at a frequency £1. We assume a tight harmonic con- 
finement along the axis of rotation such that the axial 
ground state energy far exceeds any other transverse en- 
ergy scale, yielding a quasi-2D system. The Hamiltonian 
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in the co- rotating frame is then H — H + H lnt , where 



and 



N 
i<j 



(1) 



Here, M is the particle mass, N the number of bosons, 
lu the radial trap frequency and g — VSwfiLj£ 2 a/£ z is the 
2D-interaction strength where variables I = \Jr\/Moj and 
t z = ^h/MuJz are the characteristic oscillator lengths 
along the radial and axial directions respectively, and a 
is the three-dimensional scattering length [31]. 

We proceed to calculate the energy spectrum of the 
rotating Bose gas by exact diagonalization of the Hamil- 
tonians H[ nt and H in blocks of definite total angular 
momentum. To do so, we choose a Fock basis of the 
form 

where 6^ creates a boson in state k (in our case the in- 
dex k represents a distinct pair of numbers (n,m), the 
principle quantum number and the projection of the an- 
gular momentum) and TVk is the occupation number of 
level k. The Hamiltonians are then written in terms of 
the Bose field operators where Ho = J ip ' (r)Hotp(r)dr 
and H int = \ J ^(r)^(r')H int tP(r')tP(r)drdr'. Expand- 
ing the field operators in terms of 2D harmonic oscillator 
basis states, fo(r) = SkfyL^ k ( r )' one obtains 

Ho = ^2 ^k£k 



and 
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where the £k = huj(2n + \m\ + 1 — mQ) are the single 
particle 2D harmonic oscillator eigenvalues, and ft ~ tl/uj 
is the dimensionless rotation frequency. 

In our Fock basis, Ho is diagonal with eigenvalues of 
E n /huj — 2J2f n i + J2f \ m i\ + N — QL, where L = 
~^2^ rrii is the total angular momentum of the system in 
units of H. Thus, in the purely non-interacting case, the 
energy spectrum consists solely of degenerate levels sep- 
arated by an energy gap 2Hlo. We will adopt the usual 
terminology where these different levels are called Lan- 
dau Levels, in reference to the similarities between the 
rotating Bose gas and the 2D electron gas subjected to 
a strong perpendicular magnetic field. The LLL approx- 
imation corresponds to enforcing that all particles have 
rrii > and m — 0, which greatly reduces the Hilbcrt 
space dimension. The interaction Hamiltonian, on the 
other hand, is not diagonal in our basis and so once it 
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FIG. 1: Energy spectrum as a function of total angular mo- 
mentum for N — 5 bosons in the Laughlin ground state, with 
three Landau Levels. The rotation rate was set to 51 = 1 and 
the interaction strength to g — 1. States with L < are not 
considered in this work since these are of higher energy than 
those with L > (due to the energy eigenvalue Doppler shift 
term —QL). 



is included the degeneracy of the Landau levels is lifted 
and each Landau level is split into a multitude of distinct 
states. The interaction strength is parameterized by the 
dimensionless coupling constant g = g/huil 2 = \f^Kaji z . 

For small number of particles it is feasible to calcu- 
late the entire spectrum, but for larger N and greater 
number of Landau Levels the Hilbert space becomes too 
large. In these cases, only the lowest-lying eigenvalues 
were calculated through the use of a Lanczos diagonal- 
ization algorithm. Even so, computational constraints 
have limited our calculations to the first two or three 
Landau Levels, which are obtained by including states 
where [n* + (\rrii\ — rrii)/ 2] < 1 or 2, respectively. 
The 5 particle spectrum with g = 1 and il = 1 is shown 
in Fig. 1, including the three lowest Landau levels. 



III. REPULSIVE INTERACTIONS 

We begin by considering a repulsive Bose gas. As the 
interaction energy is increased, mixing between the differ- 
ent Landau Levels becomes more important and at some 
critical interaction strength g max , the LLL approximation 
can no longer accurately describe the rotating gas. The 
standard criterion for the validity of the LLL approxima- 
tion is that the interaction energy should be smaller than 
the spacing between Landau Levels, or gp < 2Tilo where 
p is the particle density [1, 32-35] In unitlcss form, the 
crossover between the weakly and strongly interacting 
regime occurs when 



(3) 



This scaling does not provide an absolute estimate on 
5max, however; the coefficient of proportionality might 
be large enough that interesting quantum Hall states 
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FIG. 2: The critical rotation frequency fl c characterizing the 
transition from the L — to L — N state for 8 particles, as 
a function of interaction strength g. The regime of validity 
of the LLL approximation is found by determining at what 
value of g a difference of more than 10% appears between Q c 
found using the LLL approximation and Q c found using more 
Landau Levels. Squares, open triangles, and closed triangles 
depict the results for one (LLL approximation), two, and three 
Landau levels, respectively. This same labeling scheme will 
be used in all subsequent figures. 



could be achieved for experimentally accessible rotation 
frequencies and interaction strengths. Most important 
though, the 1/N scaling is found to be invalid for the 
most interesting ground states, as discussed below. 

For a non-rotating trap, the stationary Bose gas with 
no angular momentum has the lowest energy. However, 
as the rotation rate is increased, states with higher L 
experience a greater Doppler shift in energy than those 
with lower L and, as a result, at specific rotation frequen- 
cies f2 c ground state transitions are observed. Almost all 
stable ground states can be described by a relationship 
of the form L = a(N — b), where a and b are integers [2]. 
We will be concentrating on the L = N state (some- 
time called the single vortex state), the L = 2(N — 1) 
state, the Pfaffian (where L = (N - l) 2 /2 for odd N and 
L = N(N — 1)/ 2 for even N), and the Laughlin state for 
which L = N(N-l). 

To characterize the transition from the weakly in- 
teracting to strongly interacting regime, we consider 
changes in the critical trap rotation frequency f2 c of these 
four ground states as higher Landau levels are included in 
the eigenvalue calculations, as shown in Fig. 2. When the 
difference between fi c obtained with the LLL approxima- 
tion and that obtained with higher Landau levels in the 
spectrum become apparent, we say that the LLL approxi- 
mation is no longer applicable. Specifically, we define the 
crossover from weakly to strongly interacting regimes to 
occur when the relative error between the results of the 
LLL approximation and those with two or three Landau 



FIG. 3: Scaling of <7 max with number of particles for the L = 
N (lower set) and L — 2(N — 1) (upper set) states. The (? ma x 
values were obtained using the method described in Fig. 2. 
Open and filled triangles depict the results for two and three 
Landau levels, respectively. The straight lines are the best 
fits through the 3LL data. 

levels exceeds a threshold around 10%: 

\n c (g)-n c ( g ) LLL \ = Q1 

tt c (g) L LL 

The interaction strength for which this occurs is defined 
as (jf m ax, and depends on the number of particles N and 
the particular transition considered. The criterion of 10% 
was chosen to reflect typical experimental uncertainty; 
we explicitly considered error thresholds of 5% and 15%, 
and the results were not significantly changed. We also 
repeated the entire calculation using the gap AE between 
the ground and first excited state as the marker; however, 
this method proved to be unreliable. While g max values 
found using AE tend to be smaller than those obtained 
with the fl c criterion, there is too much scatter to obtain 
a useful relation between g max and N. This is a conse- 
quence of the small number of particles considered: as N 
is increased, the energy gap undergoes considerable fluc- 
tuations as new states are interjected into the spectrum. 
Furthermore, it is more valuable from an experimental 
perspective to consider variations in fl c since the trap ro- 
tation frequency can be directly controlled, whereas AE 
is indirectly controlled through the choice of g and f2. 

We now proceed to determine g ma x{N) for various N 
values, for each of the four ground states defined in the 
previous paragraph. ^From Fig. 3, it is apparent that, 
for the L = N state at least, 3 Landau levels (3LL's) 
or more are required to accurately describe the g ma , x {N) 
for N > 10 since the 2LL curve slightly deviates from 
the power law for these large iV's. While it would have 
been advantageous to examine the effect of even higher 
Landau levels, computational constraints have prevented 
us from doing so. Nevertheless, we are confident that 
3LL's are sufficient for our purposes since we do not ob- 
serve any significant deviation from a power law scaling 
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FIG. 4: Scaling of <? max with number of particles for the Pfaf- 
fian (lower set) and Laughlin (upper set) states. The g max 
values were obtained using the method described in Fig. 2. 
Open and filled triangles depict the results obtained by using 
two and three Landau levels, respectively. The straight lines 
are the best fits through the 2LL data for the Pfaffian state 
and 3LL data for the Laughlin state. 



for the small number of particles that we have consid- 
ered. Both the L = N and L = 2(N — 1) states give a 
scaling relation that is consistent with Eq. (3) since we 
find that g max oc A -1 . Most important, the prefactor 
for the L = 2(N — 1) transition is more than an order 
of magnitude larger than that for the single vortex state. 
The numerics indicate that this prefactor continues to in- 
crease for all the weakly correlated ground states (those 
with angular momentum L oc N), though we don't have 
enough data to give quantitative predictions. 

The results for g max (N) are markedly different for the 
Pfaffian and Laughlin states. A cursory glance at Fig. 4 
shows that in fact g max increases with particle number, 
which is a drastic departure from the scaling relation (3). 
Again, the prefactor increases as the angular momentum 
of these strongly correlated ground states (where L oc 
A 2 ) increases. The ground state with L = N = 4 is 
believed to be close to the Pfaffian [11] and is plotted with 
the other Pfaffian data, though numerically it appears to 
be mixed with the single-vortex state. Our calculations 
are restricted to A < 10 for the Pfaffian and A < 7 
for the Laughlin state, and thus it is difficult to make 
quantitative predictions based on finite-size scaling with 
such few atoms. Despite these limitations, however, the 
numerics clearly indicate that the standard criteria for 
the validity of the LLL approximation does not apply to 
highly correlated ground states. 

Now that we have established the limits on the inter- 
action strength required to ensure the validity of the LLL 
approximation (for each value of A considered), we turn 
our attention to the values of f2 c and the energy gap AE 
between the ground and first excited state as the num- 
ber of particles is increased. The central motivation is 
to use finite-size scaling to make predictions for the crit- 
ical rotation frequencies and gap sizes for values of A 
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FIG. 5: The critical frequency £l c and gap A_E are shown as 
a function of particle number N for the L = N state and for 
the L — 2(N — 1) state in the inset (which has the same axis 
labels as the main figure, suppressed for clarity). For each 
value of N the coupling constant IS (/max obtained from the 
best fit to the data shown in Fig. 3. 



approaching those relevant to future experiments. 

Using the 5 max (A) oc TV"- 1 ' 046 and TV" 101 obtained 
from the best fits to the data shown in Fig. 3, the scaling 
of il c and AE with particle number is shown in Fig. 5 
for the L = A and L = 2(N - 1) states. The results for 
the L = A state are mostly consistent with expectations 
in the sense that in the thermodynamic limit (A — > oo) 
the critical frequency does approach a constant value, 
albeit larger than the mean-field result (« 0.85 versus 
fl c = l/\/2 [36, 37]), while the gap goes to a constant 
AE w 0.2 in units of Tiio. The asymptotic behavior of 
the L — 2(N — 1) state is more interesting, with the 
large- value of fi c appearing to approach unity while 
the gap appears to close. 

The results for fi c and AE associated with the Pfaffian 
and Laughlin states are shown in Figs. 6 and 7 respec- 
tively. The results for both g — const, and the best fit 
values of g = g maK (N) are shown. The constant g is cho- 
sen to be the maximum value satisfying the LLL criterion 
for all the N in the series (which corresponds to that for 
the smallest N). The varying g values are the best fits 
to the data shown in Fig. 4. 

The data shown in Fig. 6 clearly indicate that when g is 
kept fixed the Laughlin gap approaches a constant in the 
large-A limit; this result is consistent with previous nu- 
merical studies of the neutral-atom Laughlin state on the 
surface of a sphere [7]. For g — g max (N = 2) = 3.12, the 
asymptotic gap is found to be AE(N » 1) m 0.2 in units 
of hui. Unfortunately, when keeping g fixed the critical 
rotation frequency approaches the experimentally chal- 
lenging limit O c — > 1. On the other hand, if g increases 
with N according to its maximum value g m ax(A), it ap- 
pears that the critical frequency levels off slightly before 
fi c = 1 for large N while AE is larger than in the LLL 
shown in the inset of Fig. 6. Although an exact 
value is impossible to obtain, it seems like the gap is ap- 
proximately 0.4?ia; for large N. This implies that with a 
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FIG. 6: The critical frequency O c and gap AE are shown as 
a function of particle number N for the Laughlin state. The 
coupling constant g — 3.12 is the same for all N, the value 
corresponding to the largest fixed g max value consistent with 
the LLL approximation for this state. The inset shows the 
variation of f2 c and AE where <? max oc TV 0,91 obtained from 
the best fits to the data shown in Fig. 4. 



FIG. 7: The critical frequency Q c and gap AE are shown 
as a function of particle number N for the Pfaffian. The 
coupling constant g = 1.62 is the same for all N, the value 
corresponding to the largest constant g m ax consistent with 
the LLL approximation for this state. The inset shows the 
variation of Q c and AE where (? max oc N obtained from 
the best fits to the data shown in Fig. 4. 



judicious choice for j > 1, the Laughlin state could be 
stabilized at lower rotation frequencies without violating 
the LLL condition. 

The data shown in Fig. 7 for the Pfaffian are less 
clear than those for the Laughlin state. For fixed g — 
9ma,x(N = 4) = 1.62, the critical frequency approaches 
the high rotation limit O c = 1 while the gap in fact 
appears to close for large N. On the other hand, in- 
creasing g proportional to N 1 - 6 according to the trend 
observed in Fig. 4 yields what appears to be a constant 
value for the gap and a critical frequency slightly lower 
than the centrifugal limit (though there is too much scat- 
ter in the data to make definite statements). Thus, like 
the Laughlin state, the Pfaffian can in principle be sta- 
bilized at experimentally accessible rotation frequencies 
without violating the LLL condition. 

The striking contrast between the weakly and strongly 
correlated states in the behavior of g m &x(N) is the result 
of a decreased sensitivity of energy levels to changes in 
interaction strength and number of Landau levels, as the 
total angular momentum increases. In the LLL approxi- 
mation, the energy spectrum is invariant (up to an overall 
scale factor) under changes in the interaction strength. 
The Laughlin state by definition remains a zero-energy 
eigenstate of the many-body Hamiltonian. Thus, the 
shift with g of all other levels must be directly propor- 
tional to their energies relative to that of the Laugh- 
lin state. In short, the sensitivity to interactions as a 
function of angular momentum must be proportional to 
the fi = 1 'yrast' line [32]. Numerically, we find that 
the relative ground-state energies for a given value of 
L > N decrease exponentially with angular momentum, 
E m in oc cxp(-7L), where 7 = 0.24, 0.16, and 0.08 for 
N = 4, 5, and 6, respectively. For 2 < L < N the 
relationship is linear [33, 38, 39]. Small angular momen- 
tum states have thus an inherently increased sensitivity 



to changes in the interaction strength. 

Low angular momentum states also have an increased 
sensitivity to the presence of higher Landau levels. We 
quantify this sensitivity by the resulting shift SE in 
ground-state energy. Like -E m i n , SE decreases with L, 
although the functional form is different: we find that 
SE(L) oc L x where x = -1.0, -1.2, -1.3 for 4, 5, and 
6 particles, respectively. Most relevant though is the 
scaling of SE with N for the weakly and strongly cor- 
related states. Recall that for each N, g max was ob- 
tained by setting a 10% limit on variations of Q c re- 
sulting from the encroaching presence of higher Lan- 
dau levels. Since il c characterizes the transition be- 
tween two states, any variation in £! c will directly re- 
sult from a mismatch in SE between these two states. 
We observe the L = N state to always transition from 
L = 0, and the Laughlin state from L = N(N — 2). 
After numerically evaluating SE(N) with g = 1 and 
O = 1, we find that [5E L=0 - 5E L=N ] (N) oc N 2 - 8 while 
[El=n(n-2) - ^Laughlin] (N) oc N~ - 054 . Let the differ- 
ence of SE between two states be called T. Should T = 0, 
then both state energies are equally shifted, and as a re- 
sult there is no change in the critical transition rotation 
frequency fi c . In other words, T characterizes the offset 
in il c compared to its LLL value. As N is increased, any 
increase in T can be countered by a concomitant decrease 
in the interaction strength, and vice versa, such that a 
constant higher Landau level influence is maintained. For 
the L = to L — N transition g must be decreased with 
N to keep T(N) unchanged while for the Laughlin tran- 
sition the reverse is true. The contrasting behavior of 
gmax(N) for the weakly and strongly correlated states is 
thus recovered. 

Though it is somewhat silly to extrapolate these small- 
N results to experimentally relevant values of N ~ 10 3 , 
it is nevertheless useful to approximately determine the 
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relevant experimental parameters required to stabilize 
the Laughlin and Pfaftian states. First, it is important 
that the excitation gap be larger than the temperature; 
for a typical temperature T ~ 10 nK one requires that 
AE/h > 200 Hz. Thus, stabilizing the Pfafhan state 
with an asymptotic gap of O.lhui would therefore require 
a large radial trap frequency on the order of kHz unless 
the temperature were strongly reduced. 

The Laughlin state with AE ss 0.4 is more accessible. 
This gap size corresponds to N = 7 and g ss 10 in the 
data presented in Fig. 4. Assuming that the N = 7 re- 
sult is already close to the large- TV asymptotic gap, which 
seems to be the case according to the inset of Fig. 6, 
then g = 10 corresponds to a scattering length of ap- 
proximately 8 fim w 1500aRb, assuming oj/2tt = 200Hz 
and oj z = 1kHz, and where aRb is the s-wave scattering 
length for 87 Rb [40]. Such a large scattering length can 
be obtained through the use of Feshbach resonances [41], 
but then the pseudopotential employed in the present 
calculations to describe the low-energy atomic collisions 
would have to be modified to explicitly take into account 
the presence of atomic pairs [10]. 
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FIG. 8: Spectrum of a boson gas of 5 particles with attractive 
interactions, using the LLL approximation. For this particu- 
lar case, g = —1.0 and fi = 1. 



IV. ATTRACTIVE INTERACTIONS 

We also have investigated the effect of having attrac- 
tive interactions between bosons. In the thermodynamic 
limit, it is known that gases of untrapped bosons are un- 
stable against collapse [42, 43] when g < 0. However, 
with the addition of a trapping potential, such gases can 
stably exist for a limited number of particles [44] be- 
cause of the finite zero-point energy of the trap. It has 
been suggested recently [29, 30] that for the special case 
of a rotating, attractive gas of 2D harmonically trapped 
bosons, there might not be an upper limit to the num- 
ber of particles in which a stable gas could exist. The 
claim is that the statistical pressure of the anyons com- 
prising the quantum Hall states can stabilize the cloud 
against collapse. The argument, however, has recently 
been disputed [45]. 

The main result of Ref. [29] is that the bare interaction 
is replaced by an effective constant: 



9cS = 9 T 



47T 



B1K 



Here, k is the Chern-Simons coefficient that relates the 
(statistical) magnetic field B to the particle density p 
through — kB = 4irp. In the rotating system, B = 
— 2mfi/?i and the effective interaction strength is then 



9efi = 9 T 



(4) 



and thus for a specific rotation frequency fi there exists a 
critical negative interaction strength g = g c = —Anhfl/p 
above which the gas can be stabilized by the so-called 
anyonic pressure. To consider the validity of this claim, 
we examine the two distinct cases where fi = u and fi < 
u. 

Let us begin by considering the high rotation frequency 
limit fi = uj. In this case, the energy spectrum depends 
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FIG. 9: Ground state energy for g — —0.1 showing the 
crossover to instability for N = 126 particles. 



solely on the interaction Hamiltonian (1). A sample spec- 
trum of the Bose gas when g < is shown in Fig. 8 for 
N = 5 and g = —1.0. The spectrum is identical to that 
for positive interactions, except that it is inverted about 
the Laughlin state. The ground state energy is now de- 
generate with that of the L — state, which can be 
obtained by directly evaluating Eq. (2) with all particles 
having n, m = 0: 



E L=0 = N + g 



N(N-l) 
4ir 



(5) 



It is important to emphasize that when fi = u>, no matter 
how much L is imparted to the system, the ground state 
energy will always be the same as for the L = state. 
Thus, we see that from Eq. 5 the gas is unstable (E < 0) 
when TV > 1 — (4ir/g). Accordingly, in the thermody- 
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namic limit, the attractive gas will never be stable even 
for the weakest interaction parameter: for large enough 
N the system always crosses over into an unstable regime, 
as seen in Fig. 9. This result is intuitively obvious: with 
attractive interactions, the ground state has L — 0, which 
obviously contains no anyons because it is non-rotating. 
Thus, the stabilization observed when N is smaller than 
its critical value is due solely to the zero point energy. 

If we now turn our attention to the situation where 
ft < u>, states with have energies that are no longer 

degenerate with the L = state due to the non-zero term 
(1 — Ct/u>)L in the Hq eigenvalue equation. In this case, 
as L is increased, so does its associated ground state en- 
ergy. It is therefore plausible that states with large L 
might have positive energy. According to this scheme, 
it would be preferable to have low Q since this would 
"lift" the high L states further away from the bottom 
of the spectrum — larger attractive interactions would be 
stabilized by lower rotation frequencies. However, this 
is contradictory to Eq. (4) which states that for fixed p 
(which is the case when L is kept constant) increasing 
fl increases the maximum attractive interaction strength 
g c that can be stabilized by the anyonic pressure. In 
other words, suppose that we decrease f2 while maintain- 
ing a constant amount of L. This will ensure that p re- 
mains constant and so, according to Eq. (4), \g c \ should 
increase with ft. However, at fixed L, the Hq energy 
is larger for smaller f2 meaning that the maximum at- 
tractive interaction strength for which the gas is stable 
should also be larger. This contradicts the predictions 
of Eq. (4), and thus we have one more reason to refute 
the claims brought forth in Ref. [29]. While it may be 
possible to stabilize an attractive boson gas consisting of 
a finite number particles in a high L state, the system 
will always become instable if N is allowed to increase 
without bounds. Other objections have been raised in 
the comment [45]. 



imation becomes invalid. For 'mean-field' states with 
L cx N, we obtain a g max which decreases with N. With 
this scaling, in the limit of large N both the excita- 
tion gap and critical rotation frequency for the L = N 
'single-vortex' state approach limiting values of around 
0.2huj and 0.85?iu;, respectively. In the same limit for 
the L = 2(N — 1) state, however, the gap closes and the 
critical frequency approaches the radial trap frequency. 
For the 'strongly correlated' states where L oc iV 2 , we 
find the remarkable result that the maximum interaction 
strength consistent with the LLL approximation actu- 
ally increases: g max oc N x , where x is a positive num- 
ber (= 1.6 for the Pfaffian and 0.91 for the Laughlin 
states). This in turn implies that the critical rotation 
frequencies required to stabilize these states could be no- 
ticeably smaller than the radial trap frequency. The gaps 
for the Laughlin and Pfaffian approach values of OAHuj 
and O.lhuj respectively. On other hand, if the interaction 
strength is chosen to be a constant for all N, the numerics 
show that in the large- N limit the critical frequency for 
the 'strongly correlated' states approaches the trap fre- 
quency, and that for the Laughlin state the energy gap 
approaches a non-zero value of 0.2, while for the Pfaffian 
it appears to close. 

We also have investigated the regime of negative in- 
teraction strengths, and conclude that in the thermody- 
namic limit, the trapped Bose gas is unstable against 
collapse for any rotation frequency, contrary to previous 
claims. 
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